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ABSTRACT 

We present deep color-magnitude diagrams for two Subaru Suprime-Cam fields in the Virgo Stellar 
Stream(VSS)/Virgo Overdensity(VOD) and compare them to a field centred on the highest concen- 
tration of Sagittarius (Sgr) Tidal Stream stars in the leading arm, Branch A of the bifurcation. A 
prominent population of main sequence stars is detected in all three fields and can be traced as faint 
as g « 24 mag. Using theoretical isochrone fitting we derive an age of 9.ltj".? Gyr, a median abun- 
dance of [Fe/H]=-0.70±g;^ dex and a heliocentric distance of 30.9 ± 3.0 kpc for the main sequence 
of the Sgr Stream Branch A. The dominant main sequence populations in the two VSS/VOD fields 
(A « 265°, B Q « 13°) are located at a mean distance of 23.3 ± 1.6kpc and have an age ~ 8.2 Gyr 
and an abundance [Fe/H]=— 0.67^0.12 ^ ex similar to the Sgr Stream stars. These statistically robust 
parameters, derived from photometry of 260 main sequence stars, are also in good agreement with the 
age of the main population in the Sagittarius dwarf galaxy (8.0 ±1.5 Gyr). They also agree with the 
peak in the metallicity distribution of 2 — 3 Gyr old M-giants, [Fe/H]« —0.6 dex, in the Sgr leading 
arm north. We then compare the results from the VSS/VOD fields with the Sgr Tidal Stream model 
by Law & Majewski based on a triaxial Galactic halo shape that is empirically calibrated with SDSS 
Sgr A-branch and 2MASS M giant stars. We find that the most prominent feature in the CMDs, 
the main sequence population at 23 kpc, is not explained by the model. Instead the model predicts 
in these directions a low density filamentary structure of Sgr debris stars at ~ 9 kpc and a slightly 
higher concentration of Sgr stars spread over a heliocentric distance range of 42 — 53 kpc. At best 
there is only marginal evidence for the presence of these populations in our data. Our findings then 
suggest that while there are probably some Sgr debris stars present, the dominant stellar population 
in the Virgo Overdensity originates from a different halo structure that has almost identical age and 
metallicity as some sections of the Sgr tidal stream. 

Subject headings: Galaxies: Individual: Name: Sagittarius, Galaxy: Abundances, Galaxies: Stellar 
Content, Galaxy: Halo, Galaxy: Structure 



1. INTRODUCTION 

The first decade of wide-field digital imaging has rev- 
olutionised the way the 3-D structure of the Milky Way 
can be mapped. The tomographic studies enabled by 
programs such as the Sloan Digital Sky Survey (SDSS, 
York et a l.l |2000|) or the Two Micron All Sky Survey 
(2MASS, Skrutskie et al. 2006[) have not only uncov- 



ered members from a new family of ultra-faint satel- 
lite galaxies but also revealed strong evidence for the 
presence of a significant amount of substructure in the 
halo of the Milky Way in the form of large scale stellar 
streams. Understanding the nature and origin of these 
satellite galaxies and their tidal debris has major im- 
plications for the validity of cosmological models (e.g. 



Kroupa et al. 2010| ). The most striking satellite galaxy 
current ly undergoing tidal disruption is the Sagittarius 
dwar f ([ibata et al.|[l9 94, 200l] iMartmez-Delgado et al. 
|2001[ Majew sEi et al.||2003| ). Tnis galaxy is l ocated ap- 
proximately 16 kpc from the Galactic center (Kunder & 



Chaboyer|2009[ ) , but the debris from the interaction with 



the Milky Way has been traced from 16-90 kp c galacto- 
centric distance across large areas of the sky | Maj ewski 
et al.||20031 INewberg et al1[2003ljBeiokurov et al.| |2006 ; 
Correnti et al.||2010p . |Layden fc Sarajedini[ Q2000p found 
evidence tor multiple epochs of star formation m the Sgr 
dwarf with the principal star formation epochs at 11, 5, 
and 0.5-3 Gyr and associated mean abundance values of 
[Fe/H]=— 1.3,— 0.7, and —0. 4, respectively. This p icture 
was subsequently refined by |Bellazzini et aL (2006) find- 
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ing that more than 80 precent of Sgr stars belong to a rel- 
atively metal-rich ([Fe/H]~ — 0.7dex) and intermediate- 
old age (8.0 ± 1.5 Gyr) po pulation called Pop A , after 
Bellazzini et al.|fll999| ) and [Monaco et al.| ( [20Q2| . 



Apart from major tidal streams, the analysis of large 
samples of tracers of halo substructure have further re- 
vealed the presence of numerous smaller stellar overden- 
sities possibly associated with unknown fainter streams 
or merger events in the halo. In particular for the re- 
gion of the Virgo constellation several detections of stel- 
lar overdensities have been reported covering the distance 
range 4-20 kpc. Most noticeable is a significant c oncen- 
tration of RR Lyrae (RRL) stars at a = 12.4 h ( Vivas 

Zinn et c 



eT"aT1[200Tl |V ivas 2002] |Vivas fe Zinn1|2003l 



2004 



ocated at ~19 



^pc f rom the Sun 



vezic et al. 12 005) 
Dutfau et al.||2006[ |Newberg et al.|[2007| |Prior et al | 
009b| ). This clump was also detected as an excess ol 
metal- poor, old F-typ e main sequence turnoff stars by 



yp 

Newberg et al. (|2002[) and as a dist i nct stellar clustering 



labeled Vir Z by jWaTSh et ah] pQ9| ). |Duffau et aL| ( |2006 ) 
showed via spectroscopy that the majority of the RRL 
and blue horizontal branch stars in the overdense region 
define a kinematically cold (a ~ 20 kms -1 ) feature cen- 
tered on a Galactic rest frame velocity Vgsr of ^100 
kms -1 , which they labeled the "Virgo Stellar Stream" 

(VSS). 

Furthermore , Juri c et al. (2008) and more recently 
|Bonaca et al . (2012) detected a density enhancement 
over 2000 sq degrees of sky toward the Virgo constella- 
tion by means of photometric parallax distance estimates 
of SDSS stars. They named this large scale feature the 
"Virgo Overdensity" (VOD). The VOD is estimated to 



2008 Vivas 



have a distance of 6 — 20 kpc ( Juric et al^ 
et al ||2008| |Keller et aL]|2Q09| > . In the following, we will 
use the term U VSS" when we are referring; to the feature 



identified kinematically and abundance- wise, and use the 
t erm "VOD" for the spatial overd ensity. 



Martmez-Delgado et al. (2004 2007) speculated that 



the Virgo stellar overdensity might be related to the com 
plex debris structure of the trailing arm of the Sagittarius 
stream, a scenario that is support ed by the Sgr Stream 
models for this re gion of the sky (|Law et al.||2005j [Fell- 



hauer et al 
has a spherical or oblate shape. 



<j i 1 1 

2006) assuming the Galactic halo potential 

Subsequent spectro- 
scopic follow-up stu dies of a sample of these RRLs by 
Prior et al. (2009b) not only found further RRLs as- 
sociated with the VSS, both in terms of velocity and 
abundance, but also revealed a population of metal-poor 
RRLs with large negative Vgsr velocities. These stars 
were taken to indicate the likely presence of a popula- 
tion of Sgr leading tidal tail stars, w hich are exp e cted to 
have such velocities in this region. |Prior et al. (2009a) 
also suggested that Sgr trailing debris may make a con- 
tribution to the population at positive Vgsr velocities, 
thou gh it was un l ikely to full y account for the VSS fea- 
ture. |Chou et al. (2007, 2010) have argued from a chemi- 



cal abundance point-ol-view tor the presence of Sgr debris 
stars in this region. 

The aim of the present paper is to provide a ro- 
bust estimate of the mean age of the stellar popula- 
tion that dominates the VSS /VOD region by means of 
deep imaging in two directions close to the 12. 4h clump. 
For that purpose we observed two Subaru Suprime- 
Cam fields covering 0.25 square degrees each, centered 



TABLE 1 

Coordinates of the Sgr Tidal Stream and the two Virgo 

FIELDS 



Field 



«2000 

(deg) 



#2000 

(deg) 



(deg) 



b 

(deg) 



A© 

(deg) (deg) 



1145+13 176.288 +13.95 250.07 +69.68 247.8 5.2 
1220-01 185.077 -01.35 286.95 +60.55 262.1 15.3 
1247-00 191.992 -00.75 301.07 +62.11 268.6 11.4 

around a 20 oo = 12 /l 20 m 18 s , £ 20 oo = -01 d 21 m 00 s and 
^2000 = 12 /l 47 m 58 s , J 2 ooo = -00 d 45 m 00 s , respectively. 
We also observed a field 19° away, in the direction of 
the highest star density of the Sagittarius B ranch A lead- 
ing ar m. All three fields were identified by | Walsh et al.j 
(2009) as having statistically significant overdensities of 
point sources. 

In Section 2 we describe the data acquisition, reduction 
and photometric calibration. The analysis of the CMDs 
is presented in Section 3. In Section 4 arguments are com- 
piled for the interpretation of our measurements in the 
context of the VOD being dominated by the Sgr Tidal 
Stream. An alternative interpretation of the CMDs nat- 
urally arises from the comparison with the predictions of 
the Sagittarius dwarf-Milky Way halo interaction model. 
This is tested in Section 5. Finally, in Section 6 we sum- 
marize and discuss our conclusions. 



DATA ACQUISITION, PHOTOMETRY AND 
CALIBRATION 

We obtained deep g, r CCD image s of the three fields 



(TableMT) using the Suprime-Cam (Miyazaki et al. 2002) 
on thesubaru Telescope during nights of 2009 Febru- 
ary 22 and 23 (PI. N. Arimoto). The Suprime-Cam 
consists of a 5 x 2 array of 2048 x 4096 CCD detec- 
tors and provides a fleld-of-view of 34' x 27' with a pixel 
scale of 0.202". The first night was photometric and 
the second night partially clear with the median seeing 
during both nights at 0.8". Each field was observed in 
a dithered series of 5 x 280 s exposures in the SDSS g 
band and 10 x 200 s in r. Data were processed using the 
pipeli ne software SDFRED dedicated to the Suprime- 
Cam flYagi et al.||2002| IQuchi et al ||2004| ). Each image 
was bias-subtracted and trimmed, fiat-fielded, distortion 
and atmospheric dispersion corrected, sky-subtracted, 
and combined in the usual manner. The astrometric cali- 
bration of each passband was based on a general zenithal 
polynomial projection derived from astrometric standard 
stars selected from online USNO catalog 



SourceExtractor (Bertin & Arnouts 1 1T996D was em- 
ployed to identify and analyse sources down to g ~ 
26 mag. The source catalogs include PSF and flux 
measurements inside two circular apertures (diameter: 
0.8"and 1.2") optimized for the seeing. We have tied our 
photometry to the SDSS photometric system by match- 
ing all the point sources from SDSS observed in our fields 
with our g-band and r-band detections. The magnitude 
zero-points with respect to the Galactic-extinction cor- 
rected SDSS photometry were obtained by calculating a 
3-<j clipped mean of the magnitude differences of the un- 
saturated SDSS stars that were detected in our Subaru 
data. The comparison of typically 100-200 SDSS stars 
in the magnitude interval 17.6 < g,r < 22.5 for each cal- 

1 http://ftp.nofs.navy.mil/data/fchpix/ 
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Fig. 1. — An example of star-galaxy separation: the aperture 
flux ratio versus total flux for all detected sources in the 1220 — 
01 field. Unsaturated stars populate well-defined vertical bands 
0.45 < go, 8—91.2 < 0.51 and 0.52 < 7*0.8— 7*1.2 < 0.59, respectively. 
Selecting objects in these two narrow regions enables a fair star- 
galaxy separation down to g ~23.7mag. Fainter objects can be 
stars or unresolved background galaxies (see Figj2|. Sources are 
color-coded with the local number density. 

ibration made these zero-points statistically robust with 
uncertainties less than 0.01 mag. 

To separate stars from background galaxies and other 
non-stellar objects, we compare the fluxes in the two dif- 
ferent apertures. The PSF geometry places stellar ob- 
jects in a well-defined, narrow region in the magnitude- 
flux ratio plane. Such a plot is shown in Fig. [I] for the 
1220—01 field. Unsaturated stars populate vertical bands 
centered at go.s — 91.2 = 0-52 and ro.s — n.2 = 0.54, re- 
spectively. Selecting only objects with values in these two 
bands ensures a fair star-galaxy separation in our images 
down to g ~ 23.5 — 24.0 mag, a limit that depends on the 
seeing. At fainter magnitudes the signal-to-noise gets 
too low and unresolved background galaxies blend into 
the region leading to an increasing level of contamina- 
tion. All color-magnitude diagrams (CMDs; (g — r)o vs. 
go) shown in our study are based on aperture photome- 
try with a 0.8" diameter aperture calibrated to the SDSS 
system and are restricted to objects that fall simultane- 
ously into both stellar flux ratio intervals as determined 
for each field individually. 

Artificial star tests were conducted to estimate the 
detection completeness. A catalogue of uniformly dis- 
tributed simulated stars was added to each image based 
on measured PSF parameters and noise properties. Pho- 
tometry of all the stars was then performed and typical 
photometric errors derived as the difference between in- 
put and computed magnitude. We calculated the frac- 
tion of simulated stars recovered at different locations 
across the image to test whether stellar crowding is 
present and how it affects the completeness of the ob- 
servations. The derived completeness curves (recovery 



rates as a function of magnitude) did not vary across the 
FoV and showed a steep gradient where the 50% source 
completeness in the photometry is approximately 0.5 mag 
brighter than the limiting magnitude of gu m ~ 26 mag. 

3. COLOR-MAGNITUDE DIAGRAMS 

Over the observed FoV of 918 sq arcmin several thou- 
sand sources satisfy the aperture flux selection. The 
three CMDs (Fig. E| show two distinct vertical plumes 
of Galactic stars reflecting the separation of the halo 
and thick disk stars at a color of (g — r)o ~ 0.3 and 
the thin dis k stars at (g — r)o ~ 1.35 (see Figj2^i or 
Chen et al. 2001[ ). Below go = 24 mag and around 
—0.2 < [g — r)o < 0.8 an increasing number of unre- 
solved background galaxies start dominating the CMDs. 
The Subaru photometry typically reaches down to a lim- 
iting magnitude of gu m ~ 26 mag but we will restrict our 
analysis to the interval 18 < g < 24. 

The CMDs were statistically decontaminated from un- 
wanted Galactic stars using model CMDs generated for 
each field with the TRILE GAL code flGirardi et al.|2005 



|Vanhollebeke et al.|2009[ ). As an example, the simulatec 
CMD of the Galactic field in the direction of 1145+13 
is shown in Fig. |2^. A comparison of the numbers of 
observed and predicted Galactic stars in the magnitude- 
color intervals 1.2 < g — r < 1.8, 20 < g < 24 finds good 
agreement (Table[2|: 



TABLE 2 

Star counts in the box 1.2 < g — r < 1.8, 20 < g < 24 



Field 



Observed TRILEGAL 



1145+13 
1220-01 
1247-00 



367 
536 
546 



395 
510 
568 



The cleaning procedure involves a process, in which 
a point source A in the CMD of 1145+13, 1220-00, or 
1247—01 is removed if there is a source B in the Galac- 
tic CMD that lies within the 3a photometry uncertainty 
ellipse of A. Source B is also removed from the control 
field catalog for the rest of the process to avoid repeated 
use. If more than one source is found in the ellipse, the 
closest to A is discarded. The decontaminated CMDs are 
shown in Fig. [3j We note that the stars in the simulated 
CMD represent a single realisation of the contamination, 
which explains why not all Galactic stars were removed. 
Cleaning efficiencies of 68 percent (1145+13), 63 per- 
cent (1220-00), and 59 percent (1247-01) were achieved 
in the "pure field" part of the CMD (1 < g - r < 2, 
20 < g < 24). For the < g — r < 1 color interval where 
the star density is higher compared to a pure field the 
cleaning process is expected to be more efficient because 
of the higher probability that an observed star matches 
with a star in the TRILEGAL CMD and thus is removed. 

3.1. Field 1145 + 13 in the Sgr Stream Leading Arm 

The tidal stream of the Sgr dwarf galaxy in the SDSS 
data visibly diverges westwards of a « 190° to gi ve a 
prominent bifurca tion known as Branches A and B (Be- 
lokurov et al.||2006b. Our field 1145 + 13 is locatedat 
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Fig. 2. — (a) Color-magnitude diagram of the Galactic field in the direction of 1145+13 simulated with the TRILEGAL 1.5 code. Different 
stellar populations are coded in color (thin disk - red; thick disk - green; halo - blue), (b-d) Reddening corrected color-magnitude diagrams 
of the three Subaru fields with representative photometric errors plotted on the left side of the CMD. (b) CMD of the field 1145 + 13 in the 
direction of the Sagittarius Stream A branch. A prominent main sequence population is visible with the turn-off at g ~ 21.0 mag. (c) CMDs 
of the field 1220—01 and (d) 1247—00, respectively. Similar to the 1145+13 field a prominent main sequence population is visible with the 
turn-off at g ~ 20.0 — 20.5 mag. Unresolved background galaxies appear as plume below g = 24 in the color interval —0.4 < g < 0.6. 



equatorial coordinates g^ooo = 176.29°, £2000 = 13.95° 
in the direction of the highest star density of th e Sagit- 
tarius Branch A as defined by |Belokurov et_al.| ( |20Q6[ ). 
In line with expectations, our CMD (Fig. 2b) exhibits 
a diagonal ridge of main sequence stars in the interval 
21 < g < 25.5, a feature that is even more pronounced 
in the statistically decontaminated version of the CMD 
(Fig. [3]). The Sgr Stream population has a measured 
MSTGcolor of (g-r) ~ 0.25. We derive the age, [Fe/H], 
and distance of the Sagittarius Stream using a maximum 



l ikelih ood method which closely follows [Frayn fc Gilmore 



2002 1 and was em ployed to analyse the CMD of Segue 3 



Fadely et al. 2011 ). The fundamental assumption is that 
the data is dominated by a single age/metallicity popula- 
tion. For the procedure, a suite of isochrones are fitted to 
a sample of stars, assigning to each a bivariate Gaussian 
probability function whose variance is set by the associ- 
ated photometric errors a g and cr( g - r )- We apply this 
analysis on all stars in our decontaminated CMD. For a 
given isochrone i we compute the likelihood 



A = Y[p({g,g - r}j\i,{g,g - r}ij,DMi), 
3 

where p is defined as: 

P({9, 9 - r}j\i, {g, g - r}ij ,DMi) = x 

27ra gj a (g _ r) . 

( \ 2 / 



(i) 



(2) 



exp 



a (g-r)j 



For each star j, {g, g— r}ij : DM^ are the magnitude, color, 
and de-reddened distance modulus values for isochrone 
i that maximize the likelihood of the entire data set 
{g,g — r}j in Equation [2] We take an approximate so- 
lution to finding the values of {g,g — r}ij and DMi by 
searching over a series of fine steps in g,g — r, and DM 



values for each isochrone. I nput isochrones are supplied 
by the Dartmouth library (Dotter et al. 2008), and lin- 
early interpolated at a step size ol 0.01 mag in the 2D 
color-magnitude space. The distance modulus is sam- 
pled over a range of 15.5 < m — M < 18.5 in steps of 
0.025 mag. To achieve the best sensitivity in the fitting 
process we used the isochrone segment between 0.5 mag 
brighter than the MSTO and 2 mag below. 

We calculate the maximum likelihood values Ci over 
a grid of isochrones, covering an age range from 5.5 
to 13.5 Gyr and metallicity range —2.5 < [Fe/H] < 
— 0.2dex. Grid steps are 0.5 Gyr in age, and 0.1 dex in 
[Fe/H]. With a grid of Ci values, we can locate the most 
likely value and compute confidence intervals by inter- 
polating between grid points. In addition to this inter- 
polation, we smooth the likelihood values over ~ 2 grid 
points in order to provide a more conservative estimate of 
parameter uncertainties. In Figure [4j we present the rela- 
tive density of likelihood values for the sample described 
above. We find the isochrone with the highest proba- 
bility has an age of 9.1 Gyr and [Fe/H] = —0.70, with 
68% and 95% confidence contours presented in the figure. 
The marginalized uncertainties (Table[3| about this most 
probable location correspond to an age of 9.1^'5Gyr, 

a metallicity of [Fe/H]= — 0.70^o'.20 ^ ex > anc ^ a distance 
modulus of DM = 17.45 ± 0.20 mag (d = 30.9 ± 3.0 kpc). 
They account for the varying photometric errors in the 
critical part 20 < g < 23 of the CMD, the differences 
in Dartmouth isochrone shapes between 5-12 Gyr, and 
isochrone model variations. Other sources of uncertain- 
ties such as the systematic errors inherent in the process 
of isochrone fitting, small errors in the photometric zero 
points and, at a more fundamental level, the underlying 
stellar physics and stellar populations used to calibrate 
the theoretical isochrones mean that the error bars are 
in general underestimates. Nevertheless, the conclusions 
we can draw from the results are unaffected. 
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TABLE 3 

Properties of the three Overdensities 



Field 



m-M 

(mag) 



(kpc) 



Age 
(Gyr) 



[Fe/H] 
(dex) 



1145+13 17.45 ±0.20 30.9 ± 3.0 9.1+^° -0.70 
1220-01 16.93 ±0.21 24.3 ± 2.5 7.9+5'.2 -0.66 
1247-00 16.73 ±0.19 22.2 ±2.1 8.5+J J -0.68 



+0.15 
0.20 
+0.20 
0.14 
+0.24 
-0.20 



Note. — Parameters inferred from fitting Dartmouth isochrones 
to the main sequence. 

The best-fitting isochrone is overplotted as a blue line 
in the left panel of Fig. [3] The associated heliocentric 
distance of the main sequence of 30.9 + 3.0 kpc is in excel- 
lent agreement with the est imated Sgr Stream dis tance 
of 29 kpc at that location ( |Belokurov et al. 2006). To 
illustrate that our results are independent of the cho- 
sen set of theoretical models we also present the Padova 
isochrone ( |Girardi et al.|2004 ) with the same age, [Fe/H], 
and distance (red curve). The slightly brighter turnoff 
and a slightly bluer lower giant branch corresponds to 
a metallicity difference of [Fe/H]^-[Fe/H]p = 0.07 dex 
and an age difference of 0.7 Gyr, well within the listed 
uncertainties derived from the Dartmouth isochrones. 
For a detailed compar ison of the Dartmouth an d Padova 
isochrones we refer to Dotter et al. (2007 2008[ ), but we 
note that the differences tor low mass stars is because of 
differences in the treatment of the equation of state and 
surface boundary conditions. These differences are addi- 
tional reasons why isochrone fitting was restricted to the 
brighter parts of the main sequence. 

There is not much information currently available in 
the literature about the age, metallicity and color of the 
leading tidal tail MSTO in this part of the Sgr Stream. 
Carlin et al.| ( |2012| presented the most recent results 
lor fields in the trailing tidal tail. Their Fig. 21 shows 
the spectroscopic MDFs of four fields along the regime 
75° < A < 130° with no significant metallicity gradient 
observed within the uncertainties. From the numbers in 
their Table 6 we calculate a mean [Fe/H] of —1.15 and a 
mean a^-pe/B] of 0.64. That la range of —1.77 <[Fe/H]< 
—0.49 is statistically in agreement with our quoted la 
range of -0.90 <[Fe/H]< -0.55. 



3.2. The VSS/VOD Fields 

The detection algorithm used by Walsh et al.| ([2009) 
to search for new ultra-faint Milky Way dwarf satellites 
flagged two positions, 1220 - 01 and 1247 - 00 (see Ta- 
ble]lt both located in the vicinity of two fields studied 
by [Martmez-Delgado et al. ( 2007[ ) and in the general di- 
rection of the Virgo Stellar Strea m /Virgo Overdensity as 
denned by excesses of RRL stars ([Vivas et al.|2001||Vivas 
™™ 1 | Vivas & Zirm|[2003l |Zmn et al.||2004[ |lvezic et al 



2002 



2005 



2002 
and 



and of F- type' main sequence stars ( |Newberg eta. 
, spanning an R.A./Decl. range of 175 u < a < 200 u 



-2.3° < S < 0.0°. The reddening corrected CMDs 
of 1220 - 01 and 1247 - 00 are shown in Fig. [2| Main 
sequences are clearly visible in the field-subtracted CMD 
from g^20 down to 24.0 mag, almost as conspicuous 
as in our Sgr Stream field. Following the same proce- 
dure as described before we established the best-fitting 
isochrones. They are shown in Fig. [3] with the corre- 
sponding maximum likelihood solutions in Fig.[4j We de- 



rive similar ages of 7.9^'.2 Gyr and 8.5^ j j Gyr and iden- 
tical metallicities of [Fe/H]=-0.66^;^ and -0.68^-jo, 
respectively. The associated heliocentric distances are 
24.3 ± 2.5 kpc and 22.2 ± 2.1 (see Table[3|. 

We estimate the strength of the main sequence pop- 
ulation in the two VOD fields from star counts within 
\g — r\ < 0.2 mag of the best-fitting isochrone and 3 mag 
down the MS from the turn-off. These numbers are listed 
in Table|4] together with the corresponding percentages 
normalised to the 249 stars found in the 1145 + 13 Sgr 
Stream field. The quoted uncertainties are the y/n stan- 
dard error. Although the VOD fields are « 13° away 
from the main ridge of the Sgr Tidal Stream we ob- 
serve only a moderate drop in MS star counts to about 
50 percent. 



Chou et al. (2010) pointed out that the lack of metal- 
rich RRL stars could mean that the entire VSS/VOD 
stellar population is metal-poor. With this point as 
a base, we investigate the possible presence of an 
old, metal-poor stellar population in the 1220—01 and 
1247—00 CMDs. For that purpose, we take the estimates 
of -1. 86 to -2.0 for the metallicity of VSS/VOD RRL 



(1 
P 



stars (Duffau et al. 
2009 ) /adopt an ai 



>se 

Dultau et a l ^ |2006 j[Newberg et al.|2007||Bell et al.|2008 
*rior et al.||2009a|) and overplot in the CMDs the cor 



2006 
oFT 



Prior et al. 



2009a An et al. 



.2 Gyr, and distance of 19 kpc 



responding Jj art mouth and Padova isochrones (Fig. [5|. 
In both CMDs, the MSTO region of the old, metal- 
poor isochrones seems to match a small group of stars at 
g « 19.7. However, the bright main sequence stars that 
would go with this feature are mostly absent for about 
0.7 — 1.0 mag below the turn-off until the two isochrones 
merge with the best-fitting isochrones. The lack of such 
bright main sequence stars in the interval 20 < g < 20.8 
that could be associated to the VSS/VOD suggests that 
the nearer population, i.e. at 19 kpc, is small at best 
when compared to the slightly more distant main se- 
quence stars at 23 kpc. 

To quantitatively assess the possibility that some of the 
stars above the MSTO (g < 20.2 in Fig. p| are from a 
smaller VOD stellar population located a few kpc closer 
along the line-of-sight we count stars around the MSTO 
with colors 0.2 < g — r < 0.5 in the decontaminated 
CMDs. The number of stars within one magnitude above 
and below the MSTO are listed in Tabled The 1220-01 
and 1247—00 fields have marginally more stars (50 ±15%) 
above the turn-off when compared to the 1145+13 field 
(34 ± 9%). If approximately 1/3 of the stars above the 
MSTO are associated to the Milky Way (i.e. leftover stars 
from the decontamination process), then the metal-poor, 
old VOD population is about 20-30% compared to the 
VOD main sequence stars in that part of the sky. 

We now discuss possible interpretations of the results 
from the CMD analysis. 

4. ARE THE VOD MAIN SEQUENCE STARS PART OF THE 
SGR TIDAL STREAM? 

Intriguingly, the statistically robust results, based on 
about 260 main sequence stars (0 < g — r < 1, g < 24) 
in the two VOD fields, are in excellent agreement with 
the age and metallicity derived for the MS stars in the 
1145+13 field selected to be in the direction of the high- 
est star density of the Sgr Branch A leading arm. This 
immediately raises the question, are the detected VOD 
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Fig. 3. — Foreground subtracted CMDs with the best-fitting Dartmouth isochrones (blue line) to the main sequence population (a) 
1145+13: has a metallicity of [Fe/H]=— 0.70 dex and an age of 9.1 Gyr. The associated heliocentric distance (modulus) is 30.9 kpc (m — M = 
17.45 mag), (b) 1220—01: has a metallicity of [Fe/H]=— 0.66 dex and an age of 7.9 Gyr. The associated heliocentric distance (modulus) is 
24.3 kpc (m — M = 16.93 mag), (c) 1247—00: has a metallicity of [Fe/H]=— 0.68 dex and an age of 8.5 Gyr. The associated heliocentric 
distance (modulus) is 22.2 kpc (m — M = 16.73 mag). The Padova isochrones with the same ages, metallicities, and distances are shown in 
red for comparison. 




-0.9 -0.7 
[Fe/H] 



1.5 -1.3 -1.1 -0.9 -0.7 -0.5 
[Fe/H] 



Fig. 4. — Smoothed maximum likelihood solutions in age-abundance space for the main sequences in the 1145 + 13 (left), 1220 — 01 
(middle) and 1247 — 00 (right) fields with contours at the 68% and 95% levels. Note that for all panels, any fit at relatively low abundances, 
e.g. —1.5, and at reasonable ages (< 11.5 Gyr) is well outside the 2cr boundary. 



main sequence stars part of the Sgr Tidal Stream? We 
find that the VOD results are in good agreement with 
the age for th e main population in t he Sagittarius dwarf 
(8.0 ±1.5 Gyr; |Bellazzini et al |2006| ) and the peak in the 
metallicity distribution function at [Fe/H] =—0.7 dex for 
the 2 — 3 Gyr old M-gia nts in the Sgr north leading arm in 
the Aq = 260° region ( |Chou et al.|[2007|) . We recognise 
that the M-giants might be biased to higher metallicities 
compared to a complete RGB sample but the Bellazzini 
et al. study is principally based on K-giants and does 
not indicate any major difference with the M-giant re- 
sults. We also note that the metallicity of the VOD MS 
stars is slightly lower when compared to the stars in the 



Sgr core where the metallicity distribution function has a 
wide spread from —1.0 to super-solar ±0 .2 and peaks at 
[Fe/H]=-0.4dex QBellazzini et al ||2008[ ). However, this 
difference can be explained wit h the reported abu ndance 
gradient along the tidal arms ( |Keller et aL"|2010| ). 

Consequently, based solely on the similarity between 
the MS population parameters derived for our VSS/VOD 
fields and those observed for the Sgr dwarf and the Sgr 
Tidal Stream stars, it would be possible to draw the con- 
clusion that the stellar population of the complex over- 
dense region around a=12.5h, £=0deg in the Virgo con- 
stellation is dominated by Sgr Tidal Stream main se- 
quence stars. However, such a conclusion ignores the ad- 
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TABLE 4 

Star Counts for the Three Main Sequences 



Field 



Counts Strength 



1145+13 249 100% 
1220-01 134 54±7% 
1247-00 132 53±7% 



Note. — Star counts within \g — r\ < 0.2 mag of the best- 
fitting isochrone and 3 mag down the MS from the turn-off. The 
third column gives the population strength relative to the main 
sequence in Sgr Stream Branch A. 



a) 22 




2.0 



Fig. 5. The same CMDs of the 1220-01 and 1247-00 fields 
The Dartmouth and Padova isochrones of an old 
jyr), metal poor ([Fe/H]=— 1.9) stellar population at 
19kpc as prescribed by the RRL stars in the VSS/VOD region are 
superimposed on the CMDs shown in Fig[3]as blue and red lines, 
respectively. In both fields, the MSTO region matches a small 
group of stars at g « 19.7. However, VSS/VOD-associated main 
sequence stars are mostly absent for about 0.7 — 1 mag below the 
turn-off. 



TABLE 5 

Number statistics of stars around the MSTO 



Field 



Above Below Percent 



1145+13 
1220-01 
1247-00 



33 
38 
33 



97 
82 
60 



34 ± 9% 
46 ± 13% 
55 ± 17% 



Note. — Stars are counted in the 0.2 < g — r < 0.5 color interval 
and within one magnitude above and below the MSTO. The fourth 
column gives the number ratios with yjn uncertainties. 

ditional information from models of the Sgr/Galaxy in- 
teraction, which provide further constraints. These will 
be discussed in the next section. We note, however, that, 
as discussed in Section 3.2, regardless of the origin of the 
dominant population, our data provide little evidence for 
the presence of a significant separate old and metal-poor 
population of stars at a distance of ~ 19 kpc. 

Finally, we like to point out that the RRLs probably 
associated with Sgr in this region do not pose a prob- 
lem. Abundances for 14 R RLs where the ass ociation 
with Sgr is unambiguous (e.g. Vivas et aL||2005 ) are low 
(([Fe/H]) = -1.76 with a = 0.22). Moreover, the six 
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Fig. 6. — Distribution of the LM10 model particles in the 180° < 
A© < 360° range with the points in the orbital latitude interval 
4.95° < Bq < 5.45° (top), 15.05° < B Q < 15.55° (centre) and 
11.15° < Bq < 11.65° (bottom) highlighted in purple. The triaxial 
model agrees well with the spatial location of the MS stars in the 
Sgr Tidal Stream field 1145 + 13 (filled circle] and the SDSS A- 
stream data (triangles) from iBelokurov et al.M2006|). For the two 
VOD fields the model predicts a small concentration of Sgr Tidal 
Stream stars at 9 kpc and in the distance intervals 46 < D < 54 kpc 
and 42 < D < 53 kpc, respectively, but no particles at 23 kpc where 
the most prominent population of main sequence stars is detected. 
From particle number counts we estimate that the Sgr Tidal Stream 
stars in the background should be more abundant in the 1247 — 00 
field than in the 1220 — 01 field by approximately a factor of 4.3. 



|Starkenburg et ah] ( |2009[ ) Group 1 stars, also likely Sgr 
objects, lie in this same part of the sky and have a similar 
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mean abundance. This abundance is consistent with the 
RRLs in the VSS/VOD region (whether they are VSS ob- 
jects or n ot). For instance, th e two large negative Vgsr 
RRLs in [Prior et ah] ( |2009b|) , which are likely Sgr ob- 
jects, have a bundances ( | Fe/Hl) = —1.57 similar to the 
Sgr RRLs of Prior et al. (2009a) and those in the Vivas 
et al. paper. In other words, we are not trying to deny 
that there is a metal-poor Sgr population, only that it is 
not dominant and that the mean abundance of Sgr RRL 
stars is a poor indicator of the mean abundance of the 
Sgr Stream, as the RRLs are likely older than the bulk 
of the Sgr stellar population. 

5. ARE THE VOD MAIN SEQUENCE STARS FROM A 
DIFFERENT HALO STRUCTURE? 

From the analysis of the deep CMDs we also gained 
information about the 3-D location of significant popu- 
lations of main sequence stars in three well-defined direc- 
tions. These geometrical anchor points allow a compar- 
ison with the currently b est Sgr Tidal Stream model by 
Law fc Majewskil fl2010[ hereafter LM10). The N-body 
simulation from Law &; Majewski is based on a complete 
all-sky view of the Sgr Stream. It is empirically cali- 
brated with SDSS A-branch and 2MASS M giant stars 
and predicts the heliocentric distances and radial veloc- 
ities of 100,000 particles stripped from the Sagittarius 
dwarf up to four orbits (approximately 3 Gyr) ago. 

In Fig. [6] we plot heliocentric distance data as a func- 
tion of orbital longitude for the simulated Sgr satellite 
debris from the best-fit triaxial model (a : b : c = 1 : 
0.99 : 0.72) within the assumptions of the simulation^] 
The locations of the MS stars in our three fields are 
shown as filled circles. To illustrate the model predic- 
tions we highlight those Sgr Stream model particles that 
fall within half a degree of each field's orbital latitude. 
The chosen width equals the observed field-of-view. The 
distance distribution of the model particles is bimodal in 
all three cases with a minor concentration at 9kpc and 
a slightly larger population between 46 < D < 54kpc 
and 42 < D < 53kpc, respectively. The model pre- 
diction agrees well with the distance we obtained from 
isochrone fitting for the MS stars in the 1145+13 field 
(top panel of Fig. pk. The good match was in fact an- 
ticipated as the LM10 Sgr Tidal Stream mod el has been 
calibrated w ith the SDSS A-branch data from Belo kurovl 
et al. ( |2006[ ) (triangles in Fig. [6k. However, no agreement 
is found between the model and the VOD main sequence 
stars (filled circle in centre and bottom panels of Fig.[6|: 
no Sgr Tidal Stream stars are expected at the distance 
of 23kpc where the maximum-likelihood fit located the 
main sequence in the 1220—01 and 1247—00 fields. If we 
adopt the LM10 triaxial model as the true representa- 
tion of the Sgr Tidal Stream in this part of the sky then 
the necessary conclusion would be that the observed MS 
stars at 23 kpc originate from a different stellar structure 
in the Milky Way halo. 

Do we find traces of Sgr Stream stars in our CMDs? 
To guide the eye we superimpose a Dartmouth isochrone 
with a metallicity of [Fe/H]=— 0.70 and an age of 8 Gyr at 
distances of 9 and 47 kpc over the 1220-01 and 1247-00 
CMDs in the panels of Fig. [7j Looking at the course 
]h leftna 



of the isochrone in both 



rand panels we notice a 



http:/ /www. astro. virginia.edu/~srm4n/Sgr/data. html 



diagonal ensemble of stars following the isochrone and 
running in parallel to the main sequence at 23 kpc from 
go = 18.5 mag down to about go = 22.5 mag. It is con- 
ceivable that these stars are indeed members of the small 
population of Sgr Stream stars at 9 kpc as predicted by 
the model. To enhance the signal we combine the CMDs 
of the two fields in Fig. [8j Now we can see a distinct 
narrow feature (20.5 < go < 22.5, 0.4 < (g - r) < 0.8) 
parallel and close to the isochrone at a distance of ap- 
proximately 10 kpc. In terms of distance, small distance 
spread and population size this filamentary structure is 
suggestive of the 9 kpc Sgr Tidal Stream predicted by the 
LM10 model. We note that this view offers an interest- 
ing and alternative interpretation of the origin of stars 
brighter than the 23 kpc MSTO as it has been discussed 
in Section 4.1. They are main sequence stars of the Sgr 
Stream at 9 kpc rather than the MSTO stars of an old, 
metal-poor population at 19 kpc. 

Moving on to the 47 kpc feature (righthand panels in 
Fig.[7]) reveals a difference between the 1220—01 and 
1247—00 CMDs. In the latter field we observe more stars 
scattered around the isochrone while the fewer stars in 
the 1220—01 CMD close to the isochrone lie systemati- 
cally above the line. This finding is again consistent with 
the picture we get from the LM10 model (see Fig. [6]) that 
predicts a population of Sgr Stream stars in the 1247—00 
direction larger by a factor of 4.3. Hence it is conceivable 
that this excess of stars are indeed representatives of the 
Sgr Stream spread over ^11 kpc (Am « 0.5) from 42- 
53 kpc along the line-of-sight. The bottom panel in Fig. [8] 
shows the scatter around the isochrone in the combined 
CMD. We attempted to fit the data with a second pop- 
ulation but the dominance of the 23 kpc population and 
the dispersed distribution of stars produced only numer- 
ically unstable solutions. Additional data obtained over 
a larger field may, however, allow further validation of 
the Sgr Tidal Stream model and thus its use in inter- 
preting our observations. In particular, Fig. 6 suggests 
that in the direction of the 1247—00 field, the Sgr Stream 
stars from the distant (42 — 53 kpc) wrap are relatively 
frequent. Observations in this vicinity with wider area 
coverage may then permit a stronger identification of the 
predicted Sgr population than is possible with current 
single field. 

6. SUMMARY AND CONCLUSIONS 

We report the detection of a prominent population of 
main sequence stars in two Subaru Suprime-Cam fields 
1220 — 01 and 1247 — 00 located in the overlapping region 
of the Virgo Stellar Stream(VSS) and the Virgo Overden- 
sity (VOD). The main sequence stars are at a heliocentric 
distance of 23.3 ± 1.6 kpc and have an age of 8.2^q J Gyr 

and a mean metallicity of [Fe/H]=— 0.67^o'i2 dex as es ~ 
timates from the best fitting isochrones. These parame- 
ters are similar to the age and metallicity we derived for 
the main stellar feature in our third Suprime-Cam field, 
1145 + 13, centered on the leading arm of the Sgr Stream, 
Branch A of the bifurcation. The Sgr Stream population 
has a MSTO color of (g — r)o ~ 0.25 and can be traced as 
faint as g « 25 mag in the CMD. From isochrone fitting 
we infer an age of 9.1±J;?Gyr and a mean metallicity 

of [Fe/H]=-0.70t£;2odex. We believe they are the first 
measurements of these quantities for the leading tidal 
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Fig. 7. — Dartmouth isochrones superimposed over the 1220 — 01 (top row) and 1247 — 00 CMD (bottom row) representing stellar 
populations at the distance of 9kpc and 47kpc, respectively. The isochrone at 9kpc (left column) has a metallicity of [Fe/H]=— 0.70 and 
an age of 8 Gyr. In both CMDs we notice groups of stars that follow the isochrone resembling the small Sgr Tidal Stream population at 
9kpc predicted by the LM10 model. This notion finds further support from the filamentous structure seen in the combined CMD (top 
panel, Fig. [51. The panels on the right show the same isochrone at a distance of 47kpc. Stars are scattered differently around the isochrone 
in the two fields, indicative for a larger population of Sgr Sream stars in the 1247 — 00 field, as is also predicted by the model. 
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Fig. 8. — Same as left and center panels in Figm but for the 
combined CMDs. (Top) we notice a distinct narrow feature below 
the isochrone at a distance of about lOkpc following closely the 
isochrone (21.2 < g < 22.5, 0.4 < g— r < 0.8). In terms of distance, 
small distance spread and population strength this feature could 
well be the 9kpc Sgr Tidal Stream predicted by the LM10 model. 
(Bottom) It is plausible that the stars scattered around the 47kpc 
isochrone are from the Sgr Tidal Stream that extends from 42 < 
D < 53kpc along the line-of-sight. 



tail in this part of the Sgr Stream. The associated he- 
liocentric distance of 30.9 ± 3.0 kpc agrees well with the 
estimated distance to the Sgr Tidal Stream in that di- 
rection of the Galactic halo as inferred from upper main 



seque nce and turn-off stars from SDSS (Belokurov et al 
2006]). 

The derived parameters for the main sequence stars in 
the two VSS/VOD fields are also a good match to the 
age of the dominant stellar population (Pop A) in the 
Sagit tarius dwarf galaxy (8.0 ± 1.5 Gyr; |Bellazzini et al. 
2006 ) as well as to the peak in the metallicity distribution 
function, [Fe/H]=-0.7dex, of the 2-3Gyr old M-giants 
in the Sgr north leading arm at the orbital longitude 
A = 260° flChou et al.|[2007| ). 

The agreement m metallicity may be taken as support- 
ing evidence that the detected VOD main seq uence stars 
are from the Sg r Tida l Stream, in line with the Martinez 
Delgado et al. ( 2007|) proposition (see also |Prior et^tl 
2009al) that the Sgr debris is a major contributor to the 
overdensity in Virgo. Although the measured metal- 
licity in the VSS/VOD fields ([Fe/H]=-0.67tg;^ dex) 
is slightly more metal poor than the stars in the Sgr 
core where the metallicity distribution function has a 
spread from —1. to super-solar +0.2 and peaks at 
[Fe/H]=-0.4dex flBellazzini et al.||2008D , the difference 
is consistent with the abundan c e grad ient along the tidal 
arms reported by Keller et al. (2010). 

However, it must be recalled that metallicity provides 
only a poor constraint as other Milky Way dwarf satel- 
lites massive enough to leave a prominent stellar feature 
in the MW halo can be expected to have similar [Fe/H] 



(Kirby et al. 2008 Meisner et al. 2012). To test the 



hypothesis that the VOD is dominated by stars from 
a different halo structure, we compared the CMDs of 
th e VSS/VOD fields with the Sgr Tidal Stream model 

by " ' - 



Law 



Majewski (2010) based on a triaxial Galactic 
halo shape that is empirically calibrated with SDSS Sgr 
A-branch and 2MASS M giant stars. In the surveyed di- 
rections the model makes precise predictions about the 
location and strength of different wraps of the Sgr Tidal 
Stream: a small number of Sgr Tidal Stream stars at a 
distance of ~9kpc and a slightly larger population dis- 
persed over the distance range 42 < D < 53 kpc (Fig. [6]). 
Overplotted isochrones on the CMDs at the prescribed 
distances (Figs. [7] and [8| indeed align with features that 
can be interpreted as the predicted Sgr Tidal Stream 
stars. The predicted number variation of the Sgr Tidal 
Stream stars between the two fields was also observed 
providing additional confidence in the Sgr Tidal Stream 
model. A more quantitative analysis of these features 
by means of isochrone fitting turned out to be unsuc- 
cessful because of the low contrast. A systematic, deep 
photometric study over a larger area will be necessary to 
collect more stars for further insight. 

The major finding from the analysis of the CMDs of 
the two selected fields in the complex VSS/VOD region is 
a prominent main sequence population at 23 kpc. Given 
these stars are absent in the currently best-fit N-body 
simulation of the tidal d isruption of the Sgr dwarf galaxy 



(Law & Majewski 



seems unlikely. T 



2010 ) an association with Sagittarius 
rese stars appear to belong to a dif- 
ferent stellar substructure in the Galactic halo. Addi- 
tionally, we find that at the VOD location this substruc- 
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ture has a similar age and metallicity as the Sgr Tidal 
Stream stars in the leading arm, Branch A of the bifur- 
cation. With a metallicity of [Fe/H]=— 0.67dex the as- 
sociated stars are also significantly more metal-rich than 
the VSS/VOD RRL stars ([Fe/H]« -1.9 dex). 
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